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ASPHERICAL LENS DESIGN AND IMAGING 


CRISTIAN E. GUTIERRE Z 
AND 

AHMAD SABRA 


Abstract. We design freeform lenses refracting an arbitrarily given incident field into a given fixed 
direction. In the near field case, we study the existence of lenses refracting a given bright object into 
a predefined image. We also analyze the existence of systems of mirrors that solve the near field and 
the far field problems for reflection. 


1. Introduction 


The goal of this paper is to show the existence of a lens that refracts a given variable unit field of 
directions, emitted from a planar source, into a predefined constant direction when the face closer 
to the source is given. We assume the lens is made of a material denser than the exterior medium. 
In order to avoid total reflection, we require the lower face of the lens, the incident field, and the 
predefined direction to satisfy the compatibility condition described in ( |3.4[ ). We use Snell's law 
of refraction in vector form, which after analysis, reduces the problem to solving a system of first 
order partial differential equations. Assuming <|3.4|) , we prove that a solution exists if and only if 
the incident field satisfies a curl condition, see < 13.13} . To minimize internal reflection, the first face 
of the lens can be chosen to be orthogonal to the incident field; this is described in Example |3.2| 

Using the far field analysis, we solve the following imaging problem in the near field case. 
Assume T is a bijective transformation between the planar source Q and a plane domain Q*. The 
goal is to find a lens -both faces to be determined- that refracts collimated radiation emitted from 
Q into Q*. The lower face of the lens satisfies a system of nonlinear partial differential equations 
| |4.3} . We find conditions on T so that a solution to the corresponding PDE exists. Once the lower 
face is found, the second face can be constructed using the results of Section |3.3| We use this 
analysis to design lenses when the image is a magnification/contraction of the planar source, see 
Sections 5.1 and 6.2 Adapting the methods from Sections [3] and |4j we solve in Section [8] the far 
field and the near field problems for reflection. In this case, the lens is replaced by a system of two 
mirrors. 

We remark that the refractive/reflective surfaces considered are not necessarily rotationally sym¬ 
metric, and the solutions found are explicit. This permits us to calculate and simulate numerically 
examples of application. 

To place our results in perspective, both from the theoretical and practical points of view, we 
mention the following. The problem of finding a convex, analytic, and symmetric lens focusing 
all rays from a point source into a point image was first solved in 1 FM87 I in 2d using a fixed point 
type argument. This result is generalized in ||Gutl3l l to 3d to construct freeform C 2 lenses that 
refract rays emitted from a point source into a constant direction or a point image; the reflection 
case is studied in l GS14al . Illumination problems with one point source are studied in I GH14 I. 
I GM131 for the refraction case and in [ CKQ991 , l lWan96 l. [KQ97I, lGS14bl for reflection. The case 
when the input rays are collimated is considered in 1GT131 .1 1GT14 1. 
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the Freeform Optics Meeting, Arlington, VA, June 2015, seet GSTS I. 
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The use and design of free form surfaces in modern optics received great attention recently due 
to applications to illumination, imaging, aerospace and biomedical engineering. For applications 
to illumination problems see 1 DLZG081 and [ QRW15 1, for aerospace engineering see |S PDvV05i |; 
and for automobile applications such as head-up displays and car driver-side mirrors design see 
IOttQS] and [ HP051J , respectively In addition, see 1FZW13 I for an extended number of related 
references. Due to recent technological advancements in ultra precision cutting, grinding, and 
polishing machines, manufacturing free form optical devices with high precision is now possible, 
see 1BMDN15I . The systems obtained enhance the performance of traditional designs and pro¬ 
vide more flexibility for designers 1 DMT13I . Moreover, they can achieve imaging tasks that are 
impossible with the symmetric designs. For example, in ICHlOl and ICH141 , the authors construct 
multiple free form mirror systems that increase the field of view of an observer and improve visual 
performance. 

The free form surfaces designed in this paper are a one parameter family given parametrically 
and therefore they may have self intersections. To obtain physically realizable surfaces when 
the input and output fields are collimated see Remark 3.4 the general case will be analyzed in 
the forthcoming paper |GS’]. In the present paper we do not consider input and output radiation 
intensities which requires other mathematical methods. For example, when a\ is optically inactive, 
this leads to solving Monge-Ampere type equations, see for example [ WMB05I Sect. 7.7], [RM02j, 
and I dcCMT l. Flowever, the results obtained here will be used in IIGS11 to design lenses, with both 
faces optically active, and where the input and output energies are taken into consideration. 

The paper is organized as follows. In Section|2j we briefly introduce the Snell's law of refraction 
in vector form. In Section |3j we precisely state and solve the far field problem. We introduce in 
Section [4] the near field imaging problem and formulate the corresponding system of PDEs ]4.3) . 
We find, in Section[5j necessary and sufficient conditions for the imaging map T to solve the system 
of PDEs when both the source and the target are in media with same refractive index. Theorem 
5.1| With the explicit formulas obtained, we next study in Section |5.11 the case when the map is 


a magnification and analyze the thickness of the lens obtained. In Sections [6] and [7j we solve the 
imaging problem when the target and the source are in media with different refractive indices. 
In this case, we find sufficient conditions for existence of local solutions, see Theorem |6.1| and 
Section [7] The same problem is solved in Section 6.2 in the plane with a different method using 
the Legendre transform as in IGutll l and when the map is a magnification. Finally, in Section [8] 
we solve the far field and the imaging problem for reflectors. 


2. Snell's law of refraction in vector form 

Suppose F is a surface in ]R 3 that separates two media I and II that are homogeneous and 
isotropic. Let V\ and v 2 be the velocities of propagation of light in the media I and II respectively. 
The index of refraction of the medium I is by definition n\ = c/v\, where c is the velocity of 
propagation of light in vacuum, and similarly the refractive index of the medium II is n 2 = c / v 2 ■ 
If a ray of lighj^jhaving direction x e S 2 , the unit sphere in ]R 3 , and traveling through the medium 
I strikes T at the point P, then this ray is refracted into the direction m e S 2 through the medium II 
according to the Snell law in vector form: 

(2.1) n\(x x v) = 112(111 x v), 

where v is the unit normal to the surface to T at P pointing towards medium II; see lLun64l 
Subsection 4.1]. It is tacitly assumed here that x ■ v > 0. We deduce from the following 

(a) the vectors x, m, v are all on the same plane (called plane of incidence); 

(b) the well known Snell's law in scalar form 

ni sin 0i = n 2 sin 0 2 , 


'Since the refraction angle depends on the frequency of the radiation, we assume light rays are monochromatic. 
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where 0 \ is the angle between x and v (the angle of incidence), 0 2 the angle between m and 
v (the angle of refraction). 

From (n\X - n 2 m) x v = 0, which means that the vector ryx - n 2 m is parallel to the normal 
vector v. If we set k = n 2 /n\, then 


( 2 . 2 ) 


X- KVl = Av, 


for some A e R. Notice that (2.21 determines A. In fact, taking dot product with v we get 


A = cos 0i - k cos 02, cos 0i = x ■ v > 0, and cos d 2 = m ■ v = yjl - k~ 2 [ 1 - (x ■ v) 2 ], so 
(2.3) A = x ■ v — k Vl - k~ 2 (1 - (x ■ v) 2 ). 


Refraction behaves differently for k < 1 and for k > 1. Indeed, by 1 GH09 . Subsection 2.1] 

(Cl) if k < 1 then for refraction to occur we need x • v > V'l - k 2 , and in \23) we have A > 0, and 
the refracted vector m is so that x ■ m > k; 

(C2) if k > 1 then refraction always occurs and we have x ■ m > 1 /k, and in l|2.3] we have A < 0. 


3. Refraction of a general field e(x) 

We begin stating the problem to be solved. We are given a C 2 surface 0 \ in ]R 3 (the bottom face 
of the lens), and a differentiable unit field e(x) = (ei(x),e 2 (x),e 3 (x)) := (e'(x),e 3 (x)), with e 3 (x) > 0 , 
that is defined for every x = (xi,x 2 ) in a plane domain Q. We want to construct a second surface 
o 2 (the top face of the lens) so that each ray emitted from (x, 0 ) with direction e(x) is refracted by 
the lens sandwiched by 01 and a 2 into a given fixed unit direction w. We assume that the medium 
below ffi has refractive index ti\, the medium above a 2 has refractive index 113 , and the material 
enclosed within the lens has refractive n 2 , such that n 2 > n\,n 2 , but n\ and n 2 are unrelated; see 
Figure [T] 

Define K\ = n 2 /ni, k 2 = n^ln 2 , and assume that the lower surface ti| of the lens is given by the 
graph of a C 2 function u. A ray emanating from each point (x, 0), x E Cl, with unit direction e(x) 
strikes 0 \ at some point P(x) = (<p(x),«(<p(x))); we assume that the map x i-» <p(x) is C 1 . 

At P(x), let Vi(x) be the unit normal to Oi pointing towards medium n 2 . We will assume that 
e(x) ■ Vi(x) > 0. Since Ki > 1, by (C2) the ray is refracted at P(x) into the lens with unit direction 
wi\(x) given from the vectorial Snell's law ( [ 2 . 2 ] by 

(3.1) e(x) - Kimi(x) = AiVi(x), 

where Ai(x) = e(x) ■ v 2 (x) - K\ 1 - k~ 2 (1 - (e(x) • v^x)) 2 ) < 0. The ray then continues traveling 
inside the lens and strikes the yet unknown surface 02 at some point Q(x). At Q(x), the incident 
ray has direction ni\ (x) and we want it to be refracted into the direction u: in medium n 3 . 

We introduce the distance 

d(x) = |P(x) - Q(x)|, 

and parametrize the unknown surface o 2 by the vector 

(3.2) /(x) := (cp(x), u(<p(x))) + d(x)m 1 (x), x = (xi,x 2 ) e Q, <p(x) = (cpi(x),<p 2 (x)). 

The purpose of this section is to construct a family of surfaces a 2 by proving the following 
existence theorem. 

Theorem 3.1. We are given a C 2 surface 0 \ parametrized by (cp(x),u((p(x))), zvith x i-> cp(x) C 1 , a C 1 unit 
field e(x) = (e'(x),e 3 (x)), and a unit direction zv. A surface ct 2 e C 2 such that the lens (cti,ct 2 ) refracts all 
rays zvith direction e(x), xeO, into zv exists if and only ifmfx) ■ zv > K 2 for all x e Cl and curle'(x) = 0. 
Moreover, ifVh(x) = e’(x)for some h, o 2 is then parametrized by f(x) = (<p(x), u(cp{x))) + d(x, C, zv) wii(x) 
zvhere m\ is given by (|3.3t and d by (|3.14). C is a constant chosen so that d > 0. 
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Figure 1. Far field model 


Proof. By (3.1) 

1 

(3.3) mAx) = — (e(x) - Ai(x)vi(x)). 

Kl 

Since k 2 < 1, from (Cl), for refraction to occur at Q(i') and to avoid total reflection, we need 
mi ■ w > k 2/ which from (3.3) is equivalent to 

(3.4) Ai(x)v\(x) ■ w < e(x) • w - K 1 K 2 , 

for all x e Q. This is a compatibility condition between e, 0 \, and zv. 

Now, applying the Snell's law at Q(x) we have that 

(3.5) mi - k 2 w = A 2 v 2 , 


where v 2 is the unit normal to <j 2 at Q(x) pointing towards medium n 2 . 

From (3.2) , to find f(x) it only remains to calculate d(x). We note, from (3.3) and (3.5) , that 
the vector e(x) - A 1 V 1 - K] k 2 zv is a multiple of the normal v 2 . Therefore this vector must be 
perpendicular to the tangent plane to cr 2 , which is equivalent to the following system of first order 
partial differential equations 

(3.6) (e(x) - Am - Ki k 2 w) ■ f Xl = 0 

(3.7) (e(x) - Am - kik 2 w) ■ f Xl = 0. 


To carry out the calculations, we will assume that d(x) is C 1 . Later on, we will show that this is 
a consequence of the assumptions on u and e. In fact, by (3.3) , we have that 

f Xl (x) = (cp(x),ll((p(x))) + 


d(x) —(e(x) - Ai vi) 

Kl 


Flence from 13.61 


0 = (e(x) - Ai vi - kik 2 w ) ■ 


{(p{x),u{(p{x))) + 


d(x) —(e(x) - Ai vi) 
Kl 


= (e(x) - Ai vi - kik 2 w\) ■ (c p(x), u(cp(x ))) + (e(x) - Ai v x ) ■ 


d(x) — (e(x) - Ai vi) 
Kl 
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- K1K2 W ■ 

= F(x) + A(x) - B(x). 


d(x) — (e(x) - A 1 V 1 ) 

Kl 


We calculate A{x). We have that (e(x) - Ai Vi) ■ (e(x) - Ai vi) = \e(x) - Ai vi | 2 = k 2 |wii | 2 = (ki) 2 , and 
therefore (e(x) - Ai Vi) • (e(x) - Ai Vi) X| = 0. So A(x) = Ki d Xl (x). 

Also, since v\{x) is the normal to 0 \ at the point {(p(x), n(tp(x))), we have that Vi • (<p(x), u(cp(x ))) x; = 
0, and hence F(x) = (e(x) - K 1 K 2 w) ■ (cp(x), u(cp(x))) . We obtain that d satisfies the following pde 

Ki d Xl - k 2 w ■ [(e(x) - Ai Vi)d(x)] Xl = - (e(x) - k x k 2 w) ■ (p(x), u(<p(x))) Xi 
which, since ip is a constant vector, can be written as 

(3.8) [(ki -k 2 w- (e(x) - Ai Vj)) d{x)] x< = - (e(x) - k 2 k 2 w) ■ (cp(x), u(cp(x))) Xi . 

Analogously, using ( |3.7) we obtain that d also satisfies the equation 

(3.9) [(ki — k 2 w (e(x) - Ai vi)) d(x)] X2 = - (e(x) - k 2 k 2 w) ■ (<p(x), u(cp(x))) X2 . 

Define 


H(x) = (ki - k 2 zv ■ ( e(x) - Ai vi)) d(x) 

G(x ) = (e(x) - k\k 2 w ) • (<p{x),u((p(x))). 

Notice that since u e C 2 , then vi e C 1 , and since e,(p £ C 1 and we assume d e C 1 , we have that H 
and G are both C 1 . We have that 

(e(x) - K\K 2 w) • (cp(x),u(cp(x))) Xi = G Xi (x) - e Xi (x) • (cp(x),u(cp(x))) 

The ray with direction e(x) strikes the surface 01 at the point (<p(x), u(<p(x))), then 

(<p(x),u(<p(x))) = p(x)e(x) + (x,0) 

for some positive function p. Since \e(x)\ = 1, then e(x) ■ e Xi (x) = 0, for i = 1,2, and therefore 

(3.10) e Xi • (cp(x), u((p(x))) = (x, 0) • e Xi (x); i = 1,2. 

Hence 


(e(x) - kik 2 w) • (<p(x),u(<p(x))) Xi = G Xi (x) - e Xi (x) • (x,0) 

= G Xi (x) - [e(x) • (x,0)] Xi + e(x) • (x, 0) Xr 


Therefore 1 3.8 1 and [3.9b become 


H Xl (x) = -G Xl (x) + [e(x) ■ (x,0)]. Tl - e(x) • (x,0) Xl 
H X2 (x) = -G X2 (x) + [e(x) ■ (x, 0 )] X2 - e(x) ■ (x,0) X2 , 


that is. 


(3.11) [H + G- e(x) ■ (x,0)] Xl (xi,x 2 ) = -ei(xi,x 2 ) 

(3.12) [H + G- e{x) ■ (x,0)] X2 (xi,x 2 ) = -e 2 (x lr x 2 ). 


Since e e C 1 , this implies that the mixed derivatives of H + G — e(x) ■ (x, 0) exist and are continuous. 
From !Rud76l Theorem 9.41], these mixed derivatives are equal. This implies that 


(3.13) 


curl (e'(x)) = (e 2 (x)) Xi - (ei(x )) X2 = 0 , 


with e'(x) = (ei(x),e 2 (x)), where we use the 2 dimensional definition of curl. Therefore there is a 
scalar function h(x\,x 2 ) such that 

Vh(x 1/ x 2 ) = e'(x 1 ,x 2 ), 


(Q simply connected). Then 13.111 implies that 


H + G- e(x) ■ (x,0) + h(x) = ip(x 2 ). 
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and from ( 3.12) we get that i^'(x 2 ) = 0, and so 

H + G - e(x) ■ (x, 0) + h{x) = C, 

with C a constant. Hence 

C - h(x) + e(x) ■ (x,0) - [e(x) - K 1 K 2 ZV] ■ ((p(x),u((p(x))) 


(3.14) 


d(x) = d(x,C,w) := 


K\ - K2 w ■ (e(x) - Ai vi(x)) 

Notice that K\ - k 2 w ■ (e(x) - Ai Vi) = K\ - KiK 2 mi • w > K\ - K\ k 2 > 0, since n 2 > n 3 . d(x,C,w) 
represents the thickness of the lens along the direction ni\(x). We obtained then a family of lenses 
with thickness depending on the choice of the constant C. This is studied more explicitly for 


magnifying lenses in Section 5.1 


At this point we have proved that if / is given by (3.2) , with d £ C 1 , and satisfying the pdes (3.6) 
and (3.7) , then curl {e.\(x,\, x 2 ), e 2 (xi, x 2 )) = 0, and d is given by the formula (3.14) with V/z = (e\, e 2 ). 
Reciprocally, if u £ C 2 , (p,e £ C 1 , with curl (ei (x), e 2 (x)) = 0 in Q simply connected, and d is given 
by (3.14) , with V/z = (c>|, <? 2 ), then d £ C 1 and the surface cr 2 parametrized by 

1 

f(x) = (<p(x), u((p{x))) + d(x) nii(x), where vi\(x) = — (e(x) - Ai v 2 ) 

Kl 


satisfies the equations (3.6) and | 3.7) This implies that the vector ni\(x) - k 2 hj is parallel to the 
normal to f(x), and therefore the surface described by f(x) refracts the incident vector ni \ (x) into 
the direction w. 

So far we obtained that / is a C 1 parametrization of the surface cj 2 . We claim that a 2 is in fact 
twice differentiable. To prove this we will show that the unit normal vector v 2 to u 2 is C 1 . By (33) 
1 

we have that v 2 = —(mi - k 2 ic), then from (Cl) it follows that 

A 2 


0 < A 2 = |A 2 | = |mi - % 2 w\ = ijl - 2k 2 7«i ■ w + k 2 . 

It is then enough to show that m\ is C 1 . We have that 

Ai = e(x) ■ Vi(x) - ki yj 1 - k~ 2 (1 - (e(x) • vi(x)) 2 ). 

Since Ki > 1, then 1 - k“ 2 (1 - (e(x) ■ Vi(x)) 2 ) > 0. We have e £ C 1 and Vi £ C 1 (since ai £ C 2 ) then 
we get Ai(x) e C 1 . Thus, m 3 £ C 1 from (3.3) . 

This completes the proof of the theorem. □ 

Remark 3.2. In the 2d case, that is, when Q is an interval, a\ and cr 2 are curves, and e(x) = 


(ei(x), e 2 (x)) is a two dimensional unit field, then Theorem 3.1 holds without any conditions on the 
derivative of ei. 

Remark 3.3. Combining lenses of the type described and using reversibility of optical paths, one 
can construct optical systems composed of two of these lenses that refract a given incoming field 
of directions e(x) into another given field of directions g(x). 


Remark 3.4. Since the surface cr 2 in Theorem 3.1 is given parametrically, it may have self¬ 
intersections. This depends on the constant C in | 3.14) , the values of ki,k 2 , and the size of 
the Lipschitz constants for the functions e,u,Du and (p, and the domains Q and Q'. Clearly, 
if u 2 has self-intersections, then it is not physically realizable. This is analyzed in detail in the 
forthcoming paper [ GSil . However, to illustrate this issue, for simplicity we analyze here the case 


+ The curl condition is very natural because the directions of light rays s (defined as the direction of the average Poynting 
vector) are parallel to the gradient VS where S is the wave front, i.e., S is solution to the eikonal equation and curl VS = 0, 
IBW591 Section 3.1.21. 
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when e(x) = zv = 1 * 3 , that is, the incoming and outgoing rays are collimated; here we have cp(x) = x. 
In this case, one can prove the following Lipschitz estimate for the refracted direction nii(x): 

(3.15) \mi(x) - mi(y)| < M(k 1 ) \Du(x) - Du(y)\, 

for all x, y e Q and M(k 1 ) a constant depending only on jci. This estimate follows by calculation 
from the following estimate for the normal: 


K'i(x) - vi(y)| < V5 \Du(x) - Du(y)\. 

This implies a Lipschitz estimate, linear in |C|, for the distance function d(x,C,iv). In fact, we 

KiK 2 ) u{x) and g(x) = K\ - K\K 2 zv ■ m 2 {x). Since 


C + vfx) 

can write d(x,C,zv) = - , with v(x) = -(1 


Six) 

Ki(l - k 2 ) < g(x) < Ki(l + k 2 ), we have 

C (g(y) - g{x)) + v(x)g(y) - v(y)g(x) 


\d(x, C, zv) - d(y, C,w)\ = 




Kj(i - x 2 y 


gWgiy) 
(g(y) H*) - 


(|C| + \v(y)\) |g(y) - g(x)|). 


From ( |3.15| ) 

IgW - g(y )| < K 1 K 2 \m\(x) - mi(y)| < Mi(ki, k 2 ) \Du(x) - Du(y)\ < Mi(k 1 , k 2 ) L d „|x - y|, 
with L Du the Lipschitz constant of Du in Q. Also 

\v(x) - v(y)\ < (1 + Kik 2 ) I u(x) - u(y)\ < (1 + KiK 2 ) L u \x - y|, 
where L„ is the Lipschitz constant of u in Q; and |z7(y)| < (1 + KiK 2 ) maxo u. Therefore 


(3.16) 


\d(x, C, iv) -d(y,C,w)\ < A(ki,k 2 )L„ + B(ki,k 2 ) 


|C| + (1 + KiK 2 ) maxu 



With this estimate we will show that if C is appropriate, then the surface cj 2 given by f(x) = 
(x,u(x)) + d(x, C, w) ni\ (x) cannot have self-intersections, i.e., / is injective. In fact, suppose that 
there are two points x + y such that f(x) = f(y). Then 


\(x, u(x)) - (y, »(y))| < d(y, C, w) |mi(y) - mi(x)| + |d(y, C, w) ni\ (x) - d(x, C, zv) mi(x)\ 
< d(y, C, zv) |mi(y) - m\(x)\ + \d(y, C, zv) - d(x, C, ip)|. 


We have 

d(y,C,zv) < 

We then obtain from | |3.15| and (3.16j that 


|C| + (1 + Ki k 2 ) maxn m 
Ki(l - k 2 ) 


1 < 


|C| + (1 + KiK 2 ) max £ ] u 
Ki(l - k 2 ) 

M(ki) 


^Ki (1 - K 2 ) 

:= a |C| + jS. 


M(ki) L Du + A(ki, k 2 ) L u + B(ki, k 2 ) 

M(ki) 


|C| + (1 + Kj k 2 ) maxu 


Ldu 


+ B(ki, k 2 ) ] L Du |C| + (1 + KiK 2 ) max u 


Ki(l - K 2 ) 


+ B(K 1 , K 2 ) Ldu + A(Ki, k 2 ) L u 


If we choose L u and L [)u sufficiently small, then B < 1/2. This implies that |C| > —. Therefore if 

2 a 

|C| < l/2n, the surface a 2 cannot have self intersections. Let C = maxn -v(x). If we choose L Du 
small enough such that C < ^ and pick C > C with |C| < then d(x, C, zv ) > 0 and the surface o 2 
is physically realizable. 
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3.1. Example 1. The case with one point source considered in BGutl31 Section 3] is a special case 
of the problem considered above. In fact, suppose that the rays with direction e(x) for x £ Q all 

(x,0) - V 

intersect at a virtual point V = (a, b, c) below the plane containing Q. Then e(x) = 


so e'(x) = 


(a, b) 


l(*,0)-V| 


which is equal to V v |(x,0) 


l(x,0) 

V\ and therefore curl c'{x) = 0 in Q. 


V| 


, and 


3.2. Example 2 . We are given the incident unit field e(x) = (ci (x), e 2 (x), e 3 (x)) = (e'(x), e 3 (x)) with 
x e Cl, where e'(x) = Vh(x), and e 3 (x) > 0. First, we are going to construct a surface oq that 
is orthogonal to the incident rays, and consequently minimizes internal reflection. Then using 
Theorem 3.1 we find the surface cr 2 such that the lens enclosed by oq and ct 2 refracts all the field 
e(x) into a given unit direction w. The surface G\ is parametrized by the vector 

g(x) = (x, 0) + A(x) e(x) = (x + A(x) e'(x), A(x) e 3 (x)), 

with A(x) a positive function to be determined so that e(x) is normal to g(x) for each reQ, This is 
equivalent to find A(x) so that e(x) is perpendicular to the tangent plane at g(x), that is, 

((x,0) + A(x) e(x)) x . ■ e(x) = 0, 

for i = 1,2. Which is in turn equivalent to e,(x) + A Xi (x) = 0, for i = 1,2, since \e(x)\ = 1. Therefore, if 
we choose A(x) = -h(x) + C, with C an appropriate constant so that A(x) > 0 for all x £ Q, we then 
obtain 0 \ orthogonal to the field e(x). 

To apply Theorem 3.1 we need to write g(x) = ( cp(x ), u(cp(x))) with cp(x) = x + A(x) e'(x), and u 
such that n(cp(x)) = A(x) e^x). By the inverse theorem, this is possible if the Jacobian of the map 
cp is not zero. In such case, we let u(z) = A (ip _ 1 (z)) e 3 (<p _ 1 (z)). Since vi(x) = e(x), then from 1 
Ai(x) = 1 - Ki, from ( |3.3] > ni\ (x) = e(x), and therefore from J3.14[ l we obtain 

C - C + K 1 K 2 w ■ {(x, 0) + A(x) e{x)) 


Then 


d{x) = 

Ki - K 1 K 2 ZV ■ e(x) 

yields the following simplified parametrization for 02 


f(x) = (x,0) + (A(x) + d(x)) e(x) = (x,0) + 


C - C + k\A(x) + K 1 K 2 ZV ■ (x,0) 
Ki - K 1 K 2 ZV ■ e(x) 


e(x). 


3.3. Example 3. Consider the special case when e(x) ~ zu = (0,0,1), i.e., the rays entering and 
leaving the lens have unit vertical direction e^. Clearly, in this case the curl condition < 3.13) is 
satisfied, and (p{x) = x. So the ray from (x, 0) strikes the surface or at the point P(x) = (x, u(x)). The 
normal at P(x) is given by 

, , (-Dn(x), 1) 

V!(X) = =/ 

Vl + \Du(x)\ 2 

and so the condition iq (x) ■ e 3 > 0 is satisfied. 

For the application that will be described later on, we now rewrite d(x) and f(x) in terms of u 
and its gradient. We have that 


Ai = e 3 • vi - ki Jl- 7 (1 - (e 3 • vi) 2 ) 


Vl + |Du(x )| 2 


and then 
with 


If 1 ) 


I 1 + 1 2 

\Du(x)\ 2 

y l K i) 



Aivi(x) = A(x) (■ -Du(x ), 1), 
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O* 





(3.17) 


A(x) : = 


1 + \Du(x)\ 2 


If 1) 


I 1 + 1 2 

|Dw(x )| 2 

y l K iJ 



Condition 13.41 becomes 
(3.18) 


A(x) < 1 - K\K 2 - 


We conclude that if u satisfies ( |3.18) , then by 13.2l, | |3.3) , and 13.141 


(3.19) 

fix) 

(3.20) 

d(x) 

(3.21) 

m i(x) 


= ( x,u(x )) + d(x)mi(x), 

(1 - Ki k 2 )h(x) + C 
Ki - k 2 (1 - A(x)) ' 

= — (A(x) Du(x), 1 - A(x)). 
K 1 


4. Imaging Problem 


We are given a bijective transformation T between the planar source Q and a plane domain 
Q*, and a > 0. Our goal now is to find two surfaces cii and a 2 such that the lens sandwiched 
by them refracts every ray emitted from (x, 0 ) with vertical direction e 3 into the point ( Tx,a ), and 
such that rays leaving 02 have direction e 3 . Using the construction from the previous section, 02 is 
determined by cti which has the form (x, u(x)) with 11 e C 2 to be calculated. 

More precisely, the ray with direction e 3 strikes cp at P(x) = (x,u(x)) and is refracted into the 
unit direction ni\ given by ( ]3.21] >. The ray then hits tr 2 at the point /(x) = (fi(x), f 2 (x), f 3 (x)), given 
in terms of u(x) by l |3.19) . It is then finally refracted into the medium n 3 with direction e 3 into the 
point (Tx, a), see Figure [2j As before, we assume that the material enclosed within the lens has 
refractive index n 2 , such that n 2 > n\, iij, but n\ and n 3 are unrelated; and we use the notation 
Ki = n 2 /ni,K 2 = n 3 /n 2 . 


Since rays leave a 2 with direction e 3 into the point ( Tx,a ), then Tx = (/i(x),/ 2 (x)). Therefore 
from |3.19[, 13.20), 13.21 ), we get that u must satisfy the following system of partial differential 
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equations 

(4.1) 


I* - * = - ( ~ IC ‘ IC2) “ W * ^ A(x)Dn(r), 
Kj - KiK 2 (1 - A(x)) 


where A(x) is given by (3.17) . To write ( |4.1[ ) more explicitly, notice that 

(4.2) A(x) (l + ^k\ + (k 2 - l)|D«(x)| 2 j = 1 - k 2 . 

Multiplying and dividing the right hand side of ( |4. 1) by + ^Jk 2 + (k 2 - 1 )|Dz/(x)| 2 j, and using 

(4.2) , we obtain that the pde < |4.1 ) can be re written as 

(1 - Ki k 2 )u(x) + C 


(4.3) 


(k 2 - KiK 2 ) yjicl + (k 2 - l)|Du(x )| 2 + k 2 (1 - KiK 2 ) 


Du(x) = 


Sx 


where Sx = Tx - x. 

Since K\ > 1, then from l |4.2) A(x) < 0. Hence the refraction condition 1 3. 18) follows immediately 
if K\K 2 < 1 , that is, when n 2 > n 3 . 

5. Solution of the imaging problem when n \ = n 3 
In this section we prove the following theorem. 

Theorem 5.1. Assume n 3 = n 3 . Let T be a C 1 bijective map between the set Q and a plane domain Q*, 
and a > 0. Then there exists a lens refracting all rays emitted from (x,0), x e Cl, with direction e 3 into the 
point ( Tx,a ) such that rays leave the lens with direction e 3 if and only if 

k 2 - 1 

(C 2 - (k 2 - 1) |S| 2 ) curl S = — S X D|S| 2 , 

where Sx = Tx - x, and C is a negative constant with |C| > |Sx| ^k 2 - 1. In this case, the lower face of the 
lens o\ is parametrized by (x, u(x)) where u is given by (5.5| , and the top face o 2 of the lens is parametrized 
by the vector f(x) = (x,u{x)) + d(x)m\(x), with d and m 2 given by and \S.21) respectively. Both 
surfaces 0 \ and o 2 are C 2 . 

Proof. In this case KiK 2 = 1 and then by ( |3.20| ), |3.21 1 we have that 

C 


(5.1) 


d(x) = -- 


k 1 - e 3 • mi 

Since d(x) > 0 and e 3 • ni\ <1 < x\, then C < 0. The pde ( ]4.3[ ) becomes 

Du(x) Sx 

= ~C' 


(5.2) 


^k 2 + (k 2 - 1)|Dh(x)| 2 


Taking absolute values it follows that 
(5.3) 


|Sx| < 


|C| 


fA 1 


Therefore to find u solving (5.2) , the constant C, and the map T must be chosen so that condition 
(5.3) is satisfied for all xeQ. When T is a magnification, condition (5.3) is related to the thickness 


of the lens which will be analyzed in Example 5.1 

To solve the pde | |5.2) , an extra condition on S is required, condition (5.6) . In fact, taking absolute 
values and squaring both sides of (5.2) we obtain that 


|D»(x )[ 2 (c 2 - (k 2 - l) |Sx| 2 ) = k 2 |Sx| 2 . 
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From (531, C 2 - (k 2 - l) |Sx [ 2 > 0, and so 

|Dw(x)| = 


Ki|Sx| 


which yields 


a /c 2 -(k 2 -1)|Sx| 2 ' 
k 2 + (k 2 - l) \Dll(x)\ 2 = 


K 2 C 2 


C 2 - (k 2 - l) |Sx| 2 ' 

Replacing this in ( |5.2) and using that C < 0, we obtain the equivalent linear system 

— Kl Sx 


(5.4) 


Du(x) = 


^C 2 - (k 2 - l) jsxf 


■= F(x) = (F 1 (x),F 2 (x)). 


If u is a C 2 solution to ( |5.4) , then the field F is conservative. Conversely if d Xl F j = r) X] F 2 in Q with 
Q simply connected, then 


(5.5) 


u{x) = u(x o) + I F(x) ■ dr. 


1‘ 


is a solution to ( |5.2) , where the integral is a line integral along an arbitrary path y from Xq to x in 
Cl 

The equality d X2 F] = d Xl F 2 yields a condition on S = (Si, S 2 ). In fact. 

Si <3*2 Si + S 2 d X2 S 2 


<9, 2 Si ^/c 2 - (k 2 - l) jsp + (k 2 - l) 


dx.Fi = -Kl - 


^C 2 - (k 2 - l) |Sx| 2 


-Kl 


(c 2 - (k 2 - 1) |S| 2 ) 


3/2 


Similarly 


<3x^2 = 


"Kl 


\ 3/2 


(^c 2 -(k 2 -i)|sf 
[(c 2 - (k 2 - 1 ) S 2 ) <3v 2 Si + (k 2 - 1 ) SiS 2 <3x 2 S 2 \. 


[(C 2 - (k 2 - l) S 2 ) <3 Vl S 2 + (k 2 - l) SiS 2 <3xjSi] . 


(C 2 - (k 2 - l) |S| 2 ) 

Therefore S must satisfy the following condition 

(C 2 - (k 2 - l) S 2 ) <3x2 Si + (k 2 - l) Si S 2 d X2 S 2 = (C 2 - (k 2 - l) S 2 ) <3 Xl S 2 + (k 2 - l) Si S 2 <3x, Si. 
Simplifying this expression yields 

/<9S 2 <9Si\ 2 n/c c t^ 1 dS 2 \ c 2 <3Si c2 dS 2 




dxi <3 x 2 


■ <3x? 


1 <9xi 


and if we set |S| = \j(S\X) 2 + (S 2 x) 2 , then the last condition is equivalent to 

K 2 — 1 

(5.6) (C 2 - (k 2 - 1) |S| 2 ) curl S = -±-— S x D|S| 2 , 

where X denotes the cross product in two dimensions. This completes the proof of the theorem. □ 
Notice that if Tx = ((a\ + l)xi, (a 2 + l)x 2 ), then Sx = (<q x\, a 2 x 2 ), curl S = 0 and 

S x D|S | 2 = (cq x lr a 2 x 2 ) x (2<qxi, 2rt 2 x 2 ) 

= 2ai«2*i*2(«i - a 2 ). 


So S satisfies condition (5.6 1 if and only if <q = a 2 ov <X\ =0 or a 2 = 0. In the following examples 
we analyze these cases. 
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Figure 3. Near Field magnification 


5.1. Example 1. Tx = (1 + a)x, a > 0 i.e. T is a magnification with factor 1 + a. In this case, we 
have that F in (15.41 is 


m = 


“Ki (X X 


■\Jc 2 - (kj - l) « 2 |x | 2 


and ( |5.4| > then reads 

(5.7) 

(5.8) 


u Xl (x) = 


ll X2 (x) = 


-K\ ax\ 

yjc 2 - ( 'k 2 - l) n 2 |x | 2 

-K\OiX\ 

yjc 2 - (k 2 - l) a 2 |x | 2 


Integrating with respect to xi yields 


u(x) = ——^—- Jc 2 - (k 2 - l)a 2 |.\f + (p(x 2 ). 
a I Kt - 1 


a(x 2 - l) 

Differentiating this with respect to x 2 and using ( |5.8) we obtain that <p'(x 2 ) = 0 and hence 

Kl ^C 2 - (k 2 - l) a 2 \x\ 2 + A. 


(5.9) 


u(x) = 


a(x 2 - 1 ) 


Notice that (5.91 implies that the graph of n is contained in the ellipsoid with equation (z - A) 2 + 

\2 " 


Kj \x\ 2 = 


C Kl 


a(x 2 - 1 ) 


. Figure 


shows the lens obtained when Tx = 2x,a = 6, « i = n 2 = 1, and 


n 2 = 1.52. 

We analyze now the thickness of the lens constructed. 
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(5.10) 


If Q is the closed ball B,,(0), then condition 15.31 is satisfied for x e B,,(0) when 

|C| 


ay < 




To analyze the thickness, we fix the value of u at a point, say at x = 0, and its normal vector Vo at 
that point. The refracted vector ni\ (0) at (0, »(0)) is determined by (3.21). From (5.1) 


(5.11) 


C = -d( 0) (ki - e 3 • mi(0)). 


d(0 ) represents the distance along the ray with direction i'll] (0) from the point (0, «(0)) to /(0) on 
the second surface. This means that once we prescribe the value of d( 0), that is, the thickness of 
the lens on the refracted ray with direction ni\ (0), the value of C is given by (5.11) . Notice that the 
value of the constant C depends also on the value of the normal Vo- 

1 

Since ri] < n 2 , by (C2) applied to a\, we have e 3 ■ m 1 > —. Therefore, 

k 1 

1 

Ki - 1 < x\ - e 3 • viAx) < Ki -. 

Kl 

For x = 0, since d( 0) > 0, from (5. TTJ we then obtain the following bounds for the constant C: 
(5.12) d( 0) (ki - 1) < |C| < d( 0) (id - . 

Notice that from <» C = - d(x ) (ki - e 3 • Hi] (x)) for all x e Q. Then from (5.12) , we obtain the 
following bound for the thickness d(x) for each x e B,.(0) 

d(0 ) —— < d(x) < d(0 ) Kl + 1 


K\ + 1 ' ' ’ ' Ki 

Now if we assume | 5.10) , we then obtain from | 5.12) that 

V k ? - 1 


ay < d(0) 


K 1 


This means that if we choose a and y arbitrarily, to have the desired refraction job we need to 
take the thickness d( 0) sufficiently large and C satisfying (5.11) . On the other hand, if we choose 
in advance the thickness d(0 ), and then pick a, y satisfying 


ay < d(0 ) 


Kl 


a/k?~ 


then choosing C in accordance with 15.111 we obtain, in view of | 5.12) , that |5.101 holds and 
therefore the solution u is given by (5.9) with the constant A determined by the value of u( 0). 


5.2. Example 2. If Tx = ((1 + ai)xi,x 2 ), then F in (5.4) is 

-K](aX],0) 


F(x) = 


^C 2 - (k* - l) a 2 x 2 x 
and T satisfies (5.3) when \a 3 X] \ < |C|/ ^Jk 2 - 1. Then solving (5.4) yields 

(5.13) u(xi,x 2 ) = Kl — -Jc 2 - (kj - 1 j a 2 x 2 + A. 

a( k 2 - lj ’ 
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More generally, if Tx = ( h(xi),x 2 ), then Sx = (h(xi) - Xj,0) satisfies <5.6 1 . If |/z(xj) - Xi| < 
|C|/ ^Jk 2 - 1, then solving | |5.4) we obtain for xq,x e Q that 


(5.14) 


u(x lr x 2 ) = u(xl,x%)+ f — 


-Ki(h(t) -t) 


dt. 


C 2 - (K 2 - 1)(/Z(f) - f) 2 


6. Solution of the imaging problem when n \ > n 3 


In this case K\K 2 < 1. We are seeking for a solution u to l|43j so that the lens (cii, 02 ) solves the 
problem described in Section|4j where o\ is given by the graph of u and 02 is parametrized by the 
vector f(x) = (x,n(x)) + d(x)m\ where d and m\ are given in | |3.20) , and ( 3.21} respectively. This 
case have potential applications to design underwater vision devices 11IGC56I . and in immersion 
microscopy ]glyl5 [ and [ WRWOO I. 

We remark the following. From ( |4.2| > and since K 1 K 2 < 1, we deduce that the denominator in 
formula 1 3.20) is positive. Since d(x) must be positive, it follows that the numerator (1 - Ki K 2 )zi(x)+C 
must be negative. Moreover, since the lens we want to construct is above the source, then a\ is 
above the plane containing Q, i.e., u must be positive. 

We will prove in this section the following theorem. 


Theorem 6.1. Let Cbe a negative constant, x 0 e Q, with Q, Q* plane domains with G simply connected, 
T : Q —> Q* bijective and C 1 , and Sx = Tx - x. Suppose that 


( 6 . 1 ) 

and 


|Sx 0 | < 


-c ^1 

Ki (Ki - K 2 ) ' 


(6.2) curl Sx = 0 

(6.3) S x D(|Sx| 2 ) = 0, 


for all x in a neighborhood ofxg; see Remar k\6.1 
Then given 


(6.4) 


0 < <5 0 < 


1 


1 - K1K2 


-C- 


\ 


■ffT 


there exists a neighborhood U ofxg, and a unique u e C 2 (U) sofring the PDE | |4.3) in U and satisfying 
u{x 0 ) = 60 and 


(6.5) 


0 < u(x) < 


-C 

1 - K 1 K 2 ' 


for every x e U. 

Proof We rewrite (|4.3) in the following form 


(" (x)+ t4^) (ki - k2) V^ 


-1 


•^k 2 (ki - k 2 ) 2 + (Ki - k 2 ) 2 (kJ - l)|Diz(x)| 2 + Ki(l - Ki k 2 ) 


Du(x)(k 1 - K 2 ) c 2 


-1 


Ki(Ki - K2) 2 
( 1 -KiK 2 ) X 
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Let' 


Ki(Ki - K 2 f 


K 1 K 2 


v(x) = (ll(x) + x _ )(Kl - k 2 ) ^k 2 - 1, Sx 
satisfying | |6.5) , with i/(%q) = 60 is equivalent to find v satisfying 

v(x)Dv(x) 


( 6 . 6 ) 


(6.7) 


Jk\{ki - K 2) 2 + \Dv(x)\ 2 + Ki(l - K 1 K 2 ) 

-—-—(Ki - k 2 ) Jk 2 - 1 <z>(x) < 0, 

1 - K 1 K 2 > 1 


Sx. Then to find u solving ((43}, 


= Sx, 


with v(x 0 ) = (§0 + --) (ki - k 2 ) Jk 2 - 1 := z 0 . 

V I-K 1 K 2 / V 

Assume this c exists. Taking absolute values in l | 6 . 6 ) yields 
( 6 . 8 ) |Sx| < |p(x)| = -v(x), 


and taking absolute values in 16.61 again, and squaring both sides, we obtain that 

(6.9) z> 2 (x) |Dz>(x)| 2 = ^ ^k 2 (ki - k 2 ) 2 + |Dz;(x)| 2 + Ki(l - K\x 2 )j |Sx| 2 . 

Let t(x) := tJk 2 (k\ - k 2 ) 2 + |Dz>(x)| 2 , so |Dz;(x)| 2 = f 2 (x) - k 2 (k 2 - k 2 ) 2 . Replacing in ( |6.9[ ) yields 

v 2 (x)(t 2 (x) - k 2 (ki - k 2 ) 2 ) = ( f(x) + Ki(l - K 1 K 2)) 2 |Sx| 2 . 

Therefore f(x) satisfies the quadratic equation 

(6.10) (v 2 (x) - |Sx| 2 ) X 2 - [2 Ki( 1 - KiK 2 )|Sx| 2 j X - K 2 [(Ki - K 2 ) 2 V 2 (x) + (1 - KlK 2 ) 2 |Sx| 2 j = 0. 

From ( |6.8) > the constant coefficient is negative and the quadratic coefficient is positive, hence < 6.10) 
has two solutions with opposite signs. The discriminant A of | |6.10) is 

A = 4k 2 v 2 (x) [(ki - k 2 ) 2 v 2 (x) - |Sx| 2 (l - k 2 )(k 2 - 1)] > 0. 


Since t(x) > 0, we obtain 


( 6 . 11 ) 


t(x) = f(x) := 


Ki( 1 - KiK 2 )|Sx | 2 + Ki|z>(x)| yj(K 1 - k 2 ) 2 v 2 (x) - |Sx| 2 (l - k 2 )(k 2 - 1) 


v 2 (x) - |Sx | 2 

and therefore f(x) can be written as a function of v and |S|. Therefore equation ( | 6 . 6 ) can be written 
in the following way 

(6.12) Dv(x) = F(x,z>(x)), 

with 

Ki(l - KiK 2 )\Sx\ 2 + Ki|z’(x)| yjix 1 - k 2 ) 2 v 2 (x) - |Sx| 2 (l - kI)(k 2 - 1) 


Sx 

(6.13) F(xMx))= — 


v 2 (x) - |Sx| 2 


+ Ki (1 - KiK 2 ) 


-(!»)■ 


where 
(6.14) 
for \x\ < 1. 


G(x) = 


Ki(l - K!K 2 )\X\ 2 + Ki ^(Ki - K 2 ) 2 - (1 - K 2 )(K 2 - l)|x| 2 


1 - |x| 2 


+ Ki(l - KiK 2 ) 


x := <p(\x\ 2 )x, 


Set F(x, z) = G 


Sx 


= 4> 


|Sx | 2 \ Sx 


= G 


We have shown that if v solves < 6 . 61 and ||6.7), then v solves ||6.12). 



:= (F 1 (x,z),F 2 (x,z)); Sx = (Six, S 2 x). 
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Conversely, we will solve 1 6.12) and from this will obtain the solution to ( 6.61 and (6.71. To this 
purpose we use a result from llHar02l Chapter 6 , pp. 117-118]. That is, if we assume that 


(6.15) 


dF] 3F 2 . dF 2 

^-(x,z) + —{x,z)F 2 (x,z) = ^-(x,z) + -fe(X'Z)Fi(x,z) 


holds for all (x, z) in a open set O, then for each (xo, Zo) € O there exists a neighborhood V of the 
point x 0 and a unique solution v to ]6.12) in V satisfying v(x 0 ) = Zq. In our case, ( ]6. 1 5) follows from 
the assumptions | |6.2) and | |6.3) . We postpone this verification for later (we will actually prove that 
l |6.15) is equivalent to l |6.2) and ]6.3) ). 

We then show that using Hartman's result, solving ( |6.12) with an appropriate initial condition, 
furnishes the solution of ( ] 6 . 6 ) and | |6.7) . In fact, as a consequence of ( ] 6 .1 ) we can pick Zo < 0 
satisfying 

(6.16) |Sx 0 | < |z 0 | < - ——(ki - 


KiK 2 


k 2 )V^ 


1 . 


Zo 


C 


2 )M - 1 1 

; \6.7\ and | | 6 . 8 ) 


-. Then with 


- Ki K 2 

in a neighborhood 


Notice that 1 6.16) holds if and only if ( |6.4) holds for bo = — 

Om - k 2 ) 

this choice of Zo, the solution v(x) to 1 6 . 12 ) with v(%o) = Zq satisfies 
V' of Xo- It remains to show that this solution v satisfies the original equation ( | 6 . 6 ) in a neighborhood 

_ _ _ _ 

of Xo- From the definition of f in | |6.11) and ( ] 6 . 1 3) we have F(x, v(x)) = ^ (f(x) + k 2 (1 - K\ k 2 )). 

Since F(x) solves l |6.10) , we get v 2 (x)(P(x) - k 2 (ki - k 2 ) 2 ) = (F(x) + Ki(l - KiK 2 )) 2 |Sx| 2 . Therefore 
from 1 6.12), |Du(x )| 2 = P(x) - k 2 (ki - k 2 ) 2 . Since F(x) > 0, we get that F(x) = f(x) and so equation 
1 6.12) is~i | 6 . 6 ) for x e V'. 

It remains to verify | |6.15) . From | |6.16) and by continuity, there is a neighborhood O of the 
point (xo,Zo) such that |Sx| < |z| for all (x,z) e O with z < 0. For all points (x,z) e O we have by 
calculation 


dFi _ 1 / dSi dS 2 

dx 2 (X ' Z) ~ z DGl ' ' :1 -- 


^(x,z) = —DGi 
dz Z 2 


dx 2 ' dx 2 
S 

Sx 


dF 2 1 /dSi dS 2 

—(x,z) = -DG 2 • —, — 
dx 1 Z \ dx 1 dx 1 

dF 2 1 = 

"3 (x, z) = DG 2 • S, 
dz z 2 


p dSi 


where DGi and DG 2 are evaluated at —; and S, —— are evaluated at x. Moreover, from 16.141, we 

Z dXj - 

have that 


DGi(x) = ( cp (|x| 2 ) + 2x\§' (|x| 2 ) , 2xix 2 (p' (|x| 2 )), 
DG 2 (x) = (2 X!X 2 (p' (|x| 2 ), qb (|x| 2 ) + 2x 2 c// (|x| 2 )). 


(9jF 

This implies that (DGi(x)-x)G 2 (x) = (DG 2 (x)-x)Gi(x), which gives that -^-F 2 = 
we obtain the identity for (x, z)eO 


()Fr 

dz 


F 1 . In addition. 


^2(x z) - —(x z) 
dx 2 {,) dx, ' 1 


(6.17) 


|Sx | 2 


dS\ 

dx 2 


dS 2 

dx. 


2 


I Sx| 2 ) 


,(?Si 


? dS_ 


(S - (S 2 x) 2 5 f + (Six)(S 2 x) 


<?X' 


<9x 


<9S 2 <?Si 


<9x 2 <9xi 
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where 


dSi dS 2 , g 

l(x) = -“z;— = -curl S 

OX 2 oX i 


2 <?Si 


](x) = (Syr) 2 -- (S 2 x) 2 — + (Six)(S 2 x) 

0X2 OX\ 


2 dS 2 \(dS 2 dS\ \ 1 ^ 


dx 2 dx\ 


= - S x D|S| 2 . 
2 


Notice that 16.15} holds if and only if 16.17} equals zero. Now from assumptions (|6.2} and (|6 .3 } , 
and since S is a constant multiple of S, we then obtain 16.15} for x in a neighborhood of Xq. 

Conversely, we shall prove that if |Sxol < |zo| with Zo < 0 and ( 6.15} holds in a neighborhood O 
of (xq, Zo) with |Sx| < |z| in O and z < 0, then this implies that \ 6 . 2 \ and ( |6.3} hold in a neighborhood 
of Xg. From (16.14} 


< p(r ) = Ki 


(1 - K\K 2 )r + A(r) 
1 — r 


+ (1 - K1K2) , 


0 < r < 1, where we have set A(r) = y()C| - k 2 ) 2 - (1 - k 2 )(k 2 - 1)r. Differentiating with respect to 
r we obtain 


<P'( r ) = 


2ki(1 - k 2 k 2 )A (r) - k 2 (1 - Kj)( k 2 - 1)(1 + r) + 2ki(ki - k 2 ) 2 


Let r = 


(6.18) 


|Sx| 2 


2(1 - r) 2 A(r) 

2 with (x,z) e O. Since ( |6.17} equals zero, we have 

z 2 (p (r) I(x) + 2cp' (r) J(x) = 0. 


Will multiply this equation by —(1 - r) 2 A(r). Notice first that 

— (1 - r) 2 A(r)(p(r) = (1 - KiK 2 )r(l - r)A(r) + (1 - r)A(r) 2 + (1 - k-\ ic 2 )(1 - r) 2 A(r) 
Kl 


and 


= (1 - KiK 2 )(1 - r)A(r) + (1 - r)A(r) 2 


—(1 - r) 2 A(r)2cp' (r) = 2(1 - k 2 k 2 )A (r) - (1 - k^)(k 2 - 1)(1 + r) + 2(ki - k 2 ) 2 . 
k 1 


Hence 1 6.18} can be written as 

0 = z 2 ((1 - KiK 2 )(1 - r)A(r) + (1 - r)A(r) 2 ) I(x) 

+ (2(1 - kiK 2 )A (r) - (1 - k 2 )(k 2 - 1)(1 + r) + 2(ki - k 2 ) 2 ) J(x) 

= (z 2 (l - KiK 2 )(1 - r)I(x) + 2(1 - kiK 2 ) J(x)j A(r) 

+ z 2 (l - r)A(r) 2 I(x) + (-(1 - k^)(k 2 - 1)(1 + r) + 2(%i - k 2 ) 2 ) J(x). 

This implies that 

(z 2 (l - k\k 2 )(1 - r) l(x) + 2(1 - KiK 2 ) J(xj) A (r) = 

z 2 (r - l)A(r) 2 f(x) + ((1 - k 2 )(k 2 - 1)(1 + r) - 2( Kl - k 2 ) 2 ) j{x) 
Squaring both sides of this identity yields 

(z 2 (l - KiK 2 )( 1 - r) f(x) + 2(1 - k 2 k 2 ) /(x)) 2 A(r) 2 

= (z 2 (r - l)A(r) 2 I(x) + ((1 - k^)(k 2 - 1)(1 + r) - 2(k 2 - k 2 ) 2 
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that is, 

(z 2 (l - KiK 2 )(l - r)I(x ) + 2(1 - K1K2) J(x))‘ ((kj - k 2 ) 2 - (1 - k^Xk 2 - l)r) 

= (z 2 ('' - 1) ((ki - k 2 ) 2 - (1 - k 2 )(k 2 - l)r)l(x) + ((1 - k 2 )(k 2 - 1)(1 + r) - 2 (ki - k 2 ) 2 ) J(x)) 2 . 
Substituting r = |Sx| 2 /z 2 yields 

((1 - kik 2 )(z 2 - |Sx| 2 )f(x) + 2(1 - kik 2 ) J(x)f ((ki - k 2 ) 2 z 2 - (1 - k,)(k 2 - 1) |Sx| 2 ) ^ = 

[(|Sx| 2 - z 2 ) ((ki - k 2 ) 2 z 2 - (1 - k%)(k 2 - 1) |Sx| 2 ) ^J(x) 

+ ((1 - k 2 )(k 2 - 1) (z 2 + |Sx| 2 ) - 2 (ki - k 2 ) 2 z 2 ) ^J(x) 

Multiplying by z 4 yields 


(6.19) 

z 2 ((1 - kik 2 )(z 2 - |Sx| 2 )I(x) + 2(1 - kik 2 ) /(x)) ((ki - k 2 ) 2 z 2 - (1 - Kj)(k 2 - 1) |Sx| 2 ) = 

[(|Sx| 2 - z 2 ) ((ki - k 2 ) 2 z 2 - (1 - k 2 )(k 2 - 1) |Sx| 2 )f(x) 

+ ((1 - k 2 )(k 2 - 1) (z 2 + |Sx[ 2 ) - 2(ki - k 2 ) 2 z 2 

Both sides of the last identity are polynomials of 8th degree in z with coefficients that depend on 
x. Comparing the coefficients of z 8 in 1 )6.19) , we get that 

(1 - KiK 2 ) 2 (Ki - k 2 ) 2 I 2 (x) = (Ki - k 2 ) 4 I 2 (x). 

Since (1 - KiK 2 ) 2 - (k 2 - k 2 ) 2 = (1 - k 2 )(1 - k|) i 4 0, it follows that f(x) = 0 in a neighborhood of Xo 
Thus replacing J(x) = 0 in ( |6 .1 9) yields 

(6.20) (2z (1 - Ki k 2 ) /(x)) 2 ((Kr - k 2 ) 2 z 2 - (1 - k 2 )(k 2 - l)|Sx| 2 ) = 

[((1 - k 2 )(k 2 - 1) (z 2 + |Sx| 2 ) - 2 (jci - k 2 ) 2 z 2 ) /(x)] 2 . 


Now both sides of 16.201 are polynomials of 4th degree in z, and comparing the coefficients of z 4 , 
we get that 


(6.21) [2(1 - KrK 2 )(Ki - k 2 )] 2 J 2 (x) = [(1 - k 2 )(k 2 - 1) - 2(ki - k 2 ) 2 ] 2 / 2 (x). 

The coefficients of / 2 (x) in ( |6.21) are different since 

[2(1 - KiK 2 )(Ki - K 2 )] 2 - [(1 - K 2 )(K 2 - 1) - 2(K1 - K 2 ) 2 ] 2 = -(1 - K 2 ) 2 (1 - K 2 ) 2 ^ 0. 

Therefore J(x) = 0 in a neighborhood of Xq. Since S is a constant multiple of S, we obtain l |6.2) and 
| |6.3) . This completes the proof of the theorem. 

□ 


6.1. Remarks on conditions | |6.2) and (6.3 1 . Condition i |6.2) is equivalent to the existence of 
function w such that S = (w Xl , w Xl ). Therefore | |6.3) is equivalent to the following quasi linear 
equation 


( 6 . 22 ) 


W XIX2 ((w Xl ) 2 - ( W X1 ) 2 ) + w Xl w Xl ( w X2X2 - w XlXl ) = 0 . 


Using Cauchy-Kowalevski's theorem, the Cauchy problem problem for this equation can be solved 
for a large class of initial data, | Joh82[ Chapter 2]. In fact, if we are given two analytic curves 
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y(s) = and T(s) = ((p(s), i p(s)) satisfying the non characteristic condition 

0 


det 


/' S' 

0 /' g' 

-(pip (p 1 — Ip 2 (pip 


+ 0, 


and z satisfies the compatibility condition z' = (p f ' + ip g', then there exists a unique solution zv 
solving \6.22\ locally and satisfying the initial conditions zv(f(s),g(s)) = z(s), w Xl (f(s),g(s)) = (p(s) 
and zv X2 (f(s), g(s)) = li’(s). In particular, we can construct transformations Tx such that Sx = Tx - x 
satisfy ( |6.2) and ( ]6.3[ ) and map the curve y into the curve T. 

Notice that S verifies condition (|5.6|), when S satisfies (16.2b and 1631. 


6.2. Solution with a different method in dimension two for a magnification. We present in this 
section a different method to find u when T is Tx = (1 + a)x, a + 0. When the dimension is two, 
equation | |4.3) becomes 

(1 - K-] K 2 )u(x) + C 


(6.23) 


(k 2 - K 1 K 2 ) + (k 2 - l)|H'(x)| 2 + k 2 ( 1 - K 1 K 2 ) 

We use the Legendre transform to change variables in | |6.23) ; we set 

zv(£) = -u(x) + x£, u'(x) = zv'(£) = x. 

Making this change of variables in (6.23) yields the ode 


u'(x) = 


ax 

k 2 - 1 


[(1 -k 1 k 2 )(zv'(Z)£-zu(Z)) + C]£ = 


K 2 — 1 


(k 2 - Ki k 2 ) ^k 2 + (k 2 - 1)£ 2 + k 2 ( 1 - KrK 2 ) 


w'(a 


which setting r(£) = (1—Ki k 2 )£ 2 -- 


— (k 2 - KiK 2 ) ^K 2 + (K 2 - 1 )£ 2 + K 2 (1 - KiK 2 ) 


can be written 


as 


r(£)zv'(£) + (K 1 K 2 -l)£zv(Z) = -C£. 

If tv < 0, then r(£) > 0 for all 6. On the other hand, if a > 0, then there is 0 >0 such that r(±P \) = 0. 
That is, zv satisfies a linear ode of form zv'(£) + P(£)zv(£) = Q(£) with 


(6.24) 


P(S) = 


(1 - Ki k 2 )£ 


QB) = 


CP(£) 

on the intervals (-oo ,-£\), (~£i, £i), (£\, +°°). Note that Q(£) = ---— then 


0 = u/(£)+P(£)ip(£) 


C 


Hence 

(6.25) 

Differentiating ip twice we obtain 


1 - KiK 2 

w(£) - 


P(£) = (zv ■ 


C 


C 


1 - K\K 2 

A 


1 - Ki K 2 


C 


1 - KiK 2 


)■ 


i - kik 2 e fp(m‘ 


(6.26) 

We have r'(£) = 2(1 - KiK 2 )^ 


zv"(£) = Ae~f p(i)di (P(£) 2 - P'(£)). 
a(x 2 - kik 2 )£ 


, and then from 1 6.24) 


P(£) 2 -P'(£) = 


+ (Kj - 1)5 2 
(1 - k iK2 ) 2 ^ 2 (1 - KiK 2 )(r(5) - £r'(£)) 


r 2 (5) 


r 2 (5) 
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-a( 1 - K\K2 ) k 2 


2 (0(k?"1) + K -1^2 


K1K2 + (1 - K1K2) -^k 2 + (k 2 - 1)£ 2 


The term between brackets is positive since KiK 2 < 1. Since a h 0, then the function P 2 (£) - P'(c) 
is never zero for £ + ±£j and has the sign of af ici jc 2 - 1). Therefore from \6.26) , zv" is never zero 
for £ + ±£i and has the sign of A a (kik 2 - 1). We then obtain that zv' is strictly monotone and 
hence it can be inverted, i.e., £ = (zv')~ l (x). We conclude that the solution u to ( |6.23| is given by 
u(x) = x(zv')~ 1 (x) - w((h/) _1 (x)). 


7. Solution of the Imaging Problem When n \ < n 3 

By reversibility of the optical paths this case reduces to the case when KiK 2 < 1 from Section [ 6 ] 
by switching the roles of n 3 and n\. Notice that KiK 2 > 1 if and only if n 3 > n\. Let T : G —> G* 
be a bijection, then T -1 : Q* —> G. By reversibility of the optical paths a lens sandwiched between 
the surfaces o\ and a 2 that images G into T(G) = G* exists if and only if the same lens images G* 
into T _ 1 (G*) = G. In the second case this means that the refractive index next to G* is /I 3 which is 
bigger that n\ the refractive index next to G, this means we are in the case from Section[ 6 ]with n 3 
and ri\ switched. Therefore from Theorem | 6 .1 1 a sufficient condition for the local existence of the 
lens is that the components of the map T -1 verify \6.2\ and l|6.3|. 


8. Solution of the far field and the near field problems for reflection 


We briefly indicate how to extend the results obtained for reflectors. 

8.1. Snell's law of reflection. If a ray with unit direction x strikes a reflective surface T at a point 
P, then this ray is reflected with unit direction m such that 

x X v = m X v, 

where v is a unit normal at P. We obtain then that 


( 8 . 1 ) x - m = 2(x ■ v)v. 

8.2. Far Field Problem for reflection. Consider the problem analyzed in Section [3j now for 
reflection, with a\ a given C 2 reflective surface given by the graph of a function u. Rays emanate 
from (x, 0) with x e G and with unit direction e(x) as in Section[3] Then proceeding similarly as in 
Section |3j we obtain the analogue of Theorem |3.1| for reflection. That is, we obtain that a surface 
02 exists, such that rays reflected by o\, with direction rri\ (x), are then reflected by cr 2 into a fixed 
unit direction zv, if and only if e’(x) = Y//(x) for some C 2 function h, and ni\(x) ■ zv < 1. In this case, 
02 is C 2 and is parametrized by the vector f(x) = (<p(x), u(<p(x))) + d(x)m\(x) with 


( 8 . 2 ) 

(8.3) 


mi(x) = e(x) - 2(e(x) ■ vi)vi. 


d(x) = 


C - li(x) + e(x) ■ (x, 0) - (e(x) - zv) ■ (cp(x), u(<p(x))) 

1 - ZV ■ Jtti(x) 


where v 1 is the unit normal to 0 \ at (cp(x), u(cp(x)). 

Moreover, as shown in Example |3.1| using the above analysis one can deduce the far field case 
considered in lGS14a[ Section 3] where rays are assumed to be emitted from a point source. 


8.3. Near Field Imaging Problem. We obtain the following analogue of Theorem |5.1| for reflection. 

Theorem 8.1. Let T be a C 1 bijective map betzveen the set G and a plane domain G*, and a > 0. Then 
there exists a system oftzvo mirrors ( 01 , 02 ) reflecting all rays emitted from (x,0), x e G, zvith direction 
e 3 = ( 0 , 0 , 1 ) into the point (Tx, a) such that rays leave the lens zvith direction e 3 if and only if 

curl S = 0, 














ASPHERICAL LENS DESIGN AND IMAGING 


21 


where Sx = Tx - x. In this case, the lower face of the lens o 1 is parametrized by (x,u(x)) where u is given 
by ( |8.7) , and the top face 02 of the lens is parametrized by the vector f(x) = (x,n(x)) + d(x)m\(x), with d 
and mi given by (|8.4) and ||8.5|, respectively. Both surfaces o\ and 02 are C 2 . 


Proof. In this case from (J8T3J and 

(8.4) 
and 

(8.5) mfx) = (0,0,1)-2 [(0,0,1) - 


d(x) = 


C 


1 -mi(x)-e 3 ' 


(-Du(x), 1) 


(-Dm(x), 1) 


2Du(x ) 


r4 


Vl + \Du(x)\ 2 V 1 + \Du(x)\ 2 \1 + |D»(x)| 2 ' 1 + |D»(x)| 2 / 

As in Section|4j we have Tx = (fi(x),f 2 (x)), where f(x) = (x, u(x)) + d(x) ni\(x), and hence 

Tx = x + CDu(x). 


Therefore u satisfies the pde 

( 8 . 6 ) 

where Sx = Tx — x. If x, xq e Q, with Q simply connected, then one obtains 


Du(x) = 


(8.7) 


u(x) = u(x 0 ) + ^ J 


Sx • dr 


if and only if curl(Sx) = 0, where the integral is a line integral along an arbitrary path C from Xo to 
x. □ 


9. Conclusion 

We have shown an explicit method to design a lens that refracts a given variable unit field of 
directions e(x) into a prescribed direction w. For the lens to exist, we show that the field must 
satisfy a curl-zero condition and a compatibility condition with w. Using this analysis we solve 
an imaging problem consisting in finding a lens that focuses two plane images A and B, one 
into another. We find sufficient conditions on the map T : A —» B, B = T(A), that guarantee the 
existence of the lens. These conditions depend on the relationships between the refractive indices. 
Similar problems for reflectors are solved. The methods used consist in solving first order systems 
of pdes, and we illustrate the results with several examples. 
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